Gaussian to Exponential Crossover in the Attenuation of Polarization Echoes in NMR 
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An ingenious pulse sequence devised by S. Zhang , B. H. Meier , and R. R. Ernst (Phys. Rev. Lett. 
69, 2149 (1992)) reverses the time evolution ("spin diffusion") of the local polarization in a dipolar 
coupled 1 H spin system. This refocusing originates a Polarization Echo whose amplitude attenuates 
by increasing the time tn elapsed until the dynamics is reversed. Different functional attenuations are 
found for a set of dipolar coupled systems: ferrocene, (CsHs^Fe, cymantrene, (CsH5)Mn(CO)3, and 
cobaltocene, (CsHs^Co. To control a relevant variable involved in this attenuation a pulse sequence 
has been devised to progressively reduce the dipolar dynamics. Since it reduces the evolution of 
the polarization echo it is referred as REPE sequence. Two extreme behaviors were found while 
characterizing the materials: In systems with a strong source of relaxation and slow dynamics, the 
attenuation follows an exponential law (cymantrene). In systems with a strong dipolar dynamics the 
attenuation is mainly Gaussian. By the application of the REPE sequence the characteristic time of 
the Gaussian decay is increased until the presence of an underlying dissipative mechanism is revealed 
(cobaltocene). For ferrocene, however, the attenuation remains Gaussian within the experimental 
time scale. These two behaviors suggest that the many body quantum dynamics presents an extreme 
intrinsic instability which, in the presence of small perturbations, leads to the onset of irreversibility. 
This experimental conclusion is consistent with the tendencies displayed by the numerical solutions 
of model systems. 

PACS Number: 75.40.Gb, 76.60.Lz, 05.40.+j, 75.10.Jm. 



I. INTRODUCTION 

Nuclear spin dynamics, often called spin diffusion, is a 
powerful tool to analyze both the local and long range 
structure in solids. Since the spin diffusion rate depends 
on the internuclear dipolar interactions, it contains useful 
information on the spatial proximity of the nuclei and on 
the dimensionality of the interaction network || . There- 
fore, there is a great interest in new pulse sequences ap- 
plicable to a broad range of systems. This is the case of 
the pulse sequence devised by S. Zhang, B. H. Meier, and 
R. R. Ernst || (ZME) which creates a local polarization 
(LP) in a homonuclear spin system (/-spins) and detects 
its later evolution. A central point of this sequence is that 
it uses the presence of rare S-spins to label each of the 
/-spins directly bonded to them. Therefore, no spectral 
resolution of /-spins is required. This sequence has been 
successfully applied to several systems with different pur- 
poses [Q-fjJ and, recently, it has been generalized ||] to 
be used under magic angle spinning (MAS) increasing its 
potential applications. On the other hand, it allows one 
to monitor the formation of a polarization echo (PE) H 
when, by external means, the LP evolution is reversed. 
This is achieved by a change of the sign of the effective 
dipolar Hamiltonian |J. The initial localized polariza- 
tion evolves into multiple-spin order during a time tn 
and then, suitable pulses switch H — > — [2]H, retracing 
the dynamics and building up a PE at t « tp + [i] tn. 

Despite the fact that the dominant part of the Hamil- 



tonian is reversed, the PE amplitude (Mpe) attenuates 
as the time tn grows. The study of this attenuation might 
provide a new characterization tool which so far has not 
been exploited. Furthermore, it should be possible to 
use the PE attenuation to obtain additional information 
on the underlying stochastic processes that destroy the 
dipolar spin coherence and inhibit the echo formation. 
Such procedure resembles the use of the Hahn's echoes 
pfj| to detect molecular diffusion In this sense, there 
is an obvious interest in experimental approaches allow- 
ing a proper characterization of this attenuation for dif- 
ferent dipolar coupled systems. While a dependence of 
Mpe on time and the spinning rate has been already re- 
ported ||, there is still a need to extend these studies 
to different systems. As a first step to use the PE at- 
tenuation as a practical tool, we should investigate its 
functional dependence on tp and, if possible, identify its 
physical origin. In principle, there could be many pos- 
sible mechanisms contributing to the attenuation of the 
PE. A non-exhaustive list includes: non inverted terms 
in the Hamiltonian, pulse imperfections, local relaxation, 
spin motion, etc. Our idea involves the use of a set of 
systems where most of these variables remain essentially 
unchanged. Within the metallocenes family we can find 
compounds which allow us to change the dipolar interac- 
tion network and independently the magnetic nature of 
the metallic nuclei, this is, we can change the strength 
of a local source of relaxation. This would allow us to 
analyze whether the dipolar interaction plays any role 
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in the attenuation and how the presence of irreversible 
processes manifests in the decay rate. 

In this work, we used the ZME sequence to study the 
PE attenuation in a set of polycrystalline samples. In 
all these systems, the molecular structure contains one 
or two cyclopentadienil rings where only one 1 H-spin is 
labeled by the rare 13 C-spin. Our results confirm, as it 
was suggested by recent experiments 0, that there are 
systems where the decay is exponential while there are 
others where it is Gaussian. However, a combination of 
both behaviors is observed in a third class of systems. 
In order to understand this and, at the same time, to 
test the role of the dipolar interaction, we devised an ad 
hoc multiple pulse sequence, based on the ZME one, to 
progressively slow down the dipolar dynamics while the 
non-invertible interactions are kept constant. 

While other mechanisms cannot be completely ruled 
out, our experimental results are consistent with the hy- 
pothesis that the dipolar dynamics controls the attenu- 
ation. In some samples, if the dynamics is sufficiently 
reduced, a crossover from the Gaussian to the exponen- 
tial attenuation law can be observed. We present exact 
numerical solutions of the spin dynamics in a cyclopen- 
tadienil ring and in a double ring molecule, which give 
further support to our hypothesis. They illustrate how 
the spin dynamics, in the presence of either small resid- 
ual interactions or strong relaxation processes, increases 
the PE attenuation. 



II. EXPERIMENTAL METHODS 

We analyzed the PE attenuation in polycrystalline 
samples of ferrocene p2| , (CsH^^Fe, cymantrene [ p"3[ , 
(C 5 H 5 )Mn(CO)3, and cobaltocene (it), (C 5 H 5 ) 2 Co and 
in a single crystal sample of ferrocene, all with natu- 
ral 13 C isotopic abundance. At room temperature, these 
compounds crystallize in a monoclinic form with space 
group P2\/a. The cyclopentadienil rings perform fast ro- 
tations around their five-fold symmetry axis leading to in- 
homogeneous 13 C spectra, where the resonance frequency 
depends on the angle between the external magnetic field 
and the molecular symmetry axis. All the NMR measure- 
ment were performed in a Bruker MSL-300 spectrometer, 
equipped with a standard Bruker CP-MAS probe. 

The ZME sequence to reverse the LP evolution is 
schematized in Fig. la. As we mentioned above, the 
central idea is to use the rare 13 C spin (^-spin) as a lo- 
cal probe to inject magnetization to one of the abundant 
1 H-spins (/-spins) and to capture what is left after the 
/-spins have evolved. In the first part, the S-spin is po- 
larized by means of a Hartmann-Hahn cross polarization 
p5[ (CP). After that, the S polarization is kept spin- 
locked for a time t$ while the remaining proton coherence 
decays to zero. Then, A) A short CP pulse of duration td 
selectively repolarize the /i-spin directly bonded to the 



5-spin, creating the initial localized state (td is the time 
when the first maximum in a simple cross polarization 
transfer occurs jl6]]); B) The LP evolves during a time 
t\ in the rotating frame under a strong spin-lock field 
with the effective Hamiltonian 
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are time averaged due to the fast rotation of the rings. 
Here, jj is the gyromagnetic factor of the /-spin, the r- k 's 
are their internuclear distances and Ojk are the angles be- 
tween internuclear vectors and the static magnetic field. 
C) A (ir/2) x pulse tilts the polarization to the labora- 
tory frame where the /-spins evolve with — [2] Ti. z during 
a time A (tt/2)_ x pulse leads the polarization back to 
the rotating frame. During the effective Hamiltonian 
— [2] H y , produces the refocusing that builds up the PE 
at t 2 « [f] (*i + t d ) = [I] t R @. D) A short CP pulse 
transfers the polarization back to S. E) The S polar- 
ization is detected while the protons are kept irradiated. 
The PE attenuation can be monitored as a function of 
by setting different values for t%. As discussed below, 
a sequence where the evolution in the laboratory frame 
precedes the refocusing in the rotating frame (Fig. lb) 
is more appropriate to compare attenuations in different 
systems. 

In order to test the effect of the dipolar dynamics in the 
PE attenuation, we replace the B) and C) parts of the 
ZME sequence by n defocusing and refocusing periods 
of ti = tji/n and t 2 = [^] t\ respectively. Therefore, the 
transfer of single-spin order into multiple-spin order, for a 
given time tji, can be gradually reduced by increasing n. 
Hereafter we will refer to this as the Reduced Evolution 
Polarization Echo (REPE) sequence. It includes extra 
refocusing periods t m /2 w [i] td/2 to compensate the 
proton spin evolution during the CP periods. This allows 
us to have an experimental point to properly normalize 
the data. The resulting multiple pulse sequence is shown 
in Fig. lc. For n — 1 the REPE sequence differs from 
the ZME one by the extra periods t m /2. From a practical 
point of view, it is important to adjust the relative phase 
and the RF amplitude of the X,— X,Y,— Y channels |T^ ] 
since small errors produce a significant signal loss. 



III. RESULTS AND DISCUSSION 

As a first step, we measured the PE amplitude as a 
function of tji for all polycrystalline samples. The nor- 
malized experimental data are shown in Fig. 2. These 
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experimental data correspond to molecules with their 
symmetry axes approximately perpendicular to the ex- 
ternal magnetic field (they were frequency selected from 
the 13 C spectrum). We used the ZME sequence with 
ujij/2n — 44.6 kHz, tc = 2 ms, ts = 1 ms, t c i = 85 
/is, for ferrocene and td = 95 /us, for cymantrene. In the 
case of cobaltocene, the short relaxation time Tf„ w 780 
fis causes a significant 13 C signal loss if long CP or spin 
lock periods are used. Therefore, we maximized the first 
polarization transfer by setting tc = td = 85 fis and we 
chose ts = 150 /is which is long enough to allow for the 
loss of proton spin coherence and short enough to prevent 
a considerable signal loss. 

Clearly, the three samples show a different functional 
dependence of the decay of the PE on tp. In the cy- 
mantrene sample the decay is clearly exponential while 
in ferrocene it is Gaussian. The decay for cobaltocene 
is more complex. Its fitting was possible only after the 
implementation of the pulse sequence of Fig. lc as ex- 
plained below. It contains both exponential and Gaus- 
sian contributions. Taking into account that in ferrocene 
the intermolecular /-/ interactions are very important 
while in cymantrene the rings stay relatively isolated 
within the experimental time scale, we set asymptotic 
values of and 0.2 respectively. Thus, the solid lines 
are fittings with the characteristic time as the only free 
parameter. In the case of cobaltocene there are two char- 
acteristic times. 

The exponential decay can be considered as a man- 
ifestation of the presence of a strong irreversible pro- 
cess characterized by a time t^. According to the Fermi 
Golden rule, the exponential decay is the signature of a 
weak interaction with each of the states in a wide band 
of a continuous spectrum jl8| . The interaction with each 
degree of freedom is weak but the overall effect results in a 
strong irreversible process. What interactions are present 
in cymantrene and cobaltocene but not in ferrocene? In 
the cymantrene sample, the strong local relaxation could 
be due to the quadrupolar nature of the Mn nucleus in 
a low symmetry environment while in cobaltocene it is 
more probably due to the paramagnetic nature of the 
Co(II). 

The sources of Gaussian decay in ferrocene and the 
Gaussian factor in cobaltocene are less obvious. When 
we take into account the magnitude of the intermolecu- 
lar dipolar interaction in these compounds, we observe 
that it grows from cymantrene to cobaltocene. Then it 
is important to evaluate if the magnitude of the dipolar 
interaction relative to h/r^ is playing any role in deter- 
mining the type of decay. 

It should be mentioned a subtle issue underlying the 
comparison of the attenuations in different systems. 
Since there is some /-spin evolution during the cross po- 
larization time, we do not have an experimental value 
for Mpe at tp — 0. Hence, the normalization of the 
signals depends on the functional decay chosen in the fit- 



tings. While in ferrocene and cymantrene systems this 
was not particularly demanding because of the obvious 
simple functional dependences they presented, this can 
be a serious restriction in a more general situation. In 
such a case, the pulse sequence sketched in Fig. lb 
should be used. Here, the amplitude at the maximum 
MpE(t2 — t m ;t\ = 0) constitutes an experimental ref- 
erence to normalize regardless of the functional law of 
the decay. The other way to overcome this problem is 
to employ the REPE sequence sketched in Fig. lc and 
described in the experimental section. This has the ad- 
ditional advantage of controlling the complexity reached 
by the dipolar order. 

We applied the REPE sequence to a ferrocene single 
crystal. Fig. 3 shows the 13 C-NMR spectrum (with the 
sequence parameters given in the caption). Each peak 
corresponds to one of the two magnetically incquivalent 
sites in the crystal. The carrier frequency was finely 
tuned to correspond to one of them. The good signal 
to noise ratio allows us to detect small fractions of the 
initial polarization. The attenuation of the PE for n = 1 
is shown in Fig. 4. The well-defined Gaussian decay con- 
firms the results obtained for the polycrystalline sam- 
ple. The compensating periods t m /2 allow us to prop- 
erly normalize the data at tp = 0. Besides, the un- 
bounded nature of the spin network justify an asymp- 
totic value of for Mp E - The solid line represents a 
fitting with M PE (t R ) = cxp(-i(i fl /r^ b ) 2 ), obtaining 
T ,nb = ( 245 ± 5 ) When we set n = 2, 8, 16 (inset 
Fig. 4), the Gaussian decay is kept but its characteristic 
time r^b grows. As we discuss below, this is in agreement 
with our hypothesis that the time scale for the Gaussian 
attenuation is controlled by the dipolar interaction it- 
self. We found that grows as r™ 6 = a ■ n + b with 
a = (54 ± 2) fis and b = (210 ± 10) fis. 

When applying the same sequence to the cobaltocene 
sample we observe a rather different behavior. While 
the dipolar interaction network is quite similar to the 
ferrocene one, the paramagnetic Co (II) atom introduces 
a strong source of relaxation. Thus, we expect a dipolar 
dynamics similar to that of ferrocene but convoluted with 
an extra relaxation mechanism. Fig. 5 shows the PE 
attenuation in the cobaltocene sample as a function of 
tp for n = l,2,5,8,16. There is a clear crossover from a 
Gaussian decay to an exponential one. This is the main 
result of this paper. The solid lines are fittings of the 
experimental data to the equation 
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Based on the results for ferrocene we set a linear depen- 
dence of on n. This leaves only three free parameters 
for the whole data set. By increasing n, Mp E reaches an 
asymptotic exponential decay with = (640 ± 20) fis. 
This shows that the Gaussian factor can be gradually 
reduced until the underlying exponential decay becomes 
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dominant. Once this regime is reached, further incre- 
ment of n does not produce any effect. Expression (||) 
must be considered as an empirical one to account for 
the experimental results. 

At this point it is important to mention the difference 
between our results and those that could be expected 
by considering the analogies of the REPE sequence with 
the Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence 
JjT[ used to identify the contributions that attenuate the 
Hahn echo. If a single echo pulse sequence is applied, 
the amplitude of the Hahn echo at time tn is propor- 
tional to exp[— tn/T-2 — {tH/T D ) 3 ]. Here T 2 represents 
some processes not inverted by the Hahn echo and t d 
the molecular diffusion process in presence of field gradi- 
ents. When the single echo pulse sequence is replaced by 
a train of n pulses spaced apart by tjj/n, the signal at 
the time tn, is proportional to 

LM J -/T2-{ t -/r D f]) n (4) 
y n n J 

It is obvious that for large n this becomes a single expo- 
nential, being the diffusion term eliminated by the train 
of pulses. One might think that a similar argument could 
be used to explain why a train of pulses eliminates the 
Gaussian factor in (|j|). However, one should understand 
the conditions for the validity of (0) to see whether equiv- 
alent conditions apply to the REPE sequence. In order to 
have a contribution affected by a train of pulses, its orig- 
inating process must be affected. This is independent of 
the particular decay law. The condition that justifies the 
multiplication of echo amplitudes in (Q) is that at each 
intermediate echo the initial conditions are recovered, al- 
though with a lower amplitude. This requires that the 
missed amplitude must not retain any correlation, clas- 
sical or quantum, with the observed state. In the case of 
molecular diffusion, which in principle might lead to clas- 
sical correlations, these are destroyed by the 7r-pulse be- 
cause the phase acquired from the field gradient accumu- 
lates excluding the backdiffusive events from contributing 
to the observed polarization. Now let us go back to the 
REPE sequence. On the basis of our experiment we have 
two categories of processes. The ones characterized by 
Tff, are not modified by the pulse sequence and hence are 
not related to the variables affected by these pulses. With 
respect to the other, we first observe that while ampli- 
tude is modified by n it does not follow the product law. 
From this last point we learn that correlations classical 
or quantum are present. The pulse sequence is intended 
to modify dipolar spin dynamics. However, other effects 
are simultaneous: Pulse imperfections and higher order 
terms of the average Hamiltonian which are not inverted. 
A major part of the pulse imperfections is accounted for 
by the renormalization of the attenuation curves. While 
the remaining effects could, in principle, affect the PE at- 
tenuation, reasonable magnitudes included in the numer- 
ical simulations in small systems (see next section) can- 



not account, by themselves, for the observed time scale 
of the decay. An effect present in the real experiments is 
the unboundedness of the spin system. This huge space 
becomes available to the spin correlations at a rate con- 
trolled by n. Since they are entangled with non inverted 
interactions, they could make it available to irreversible 
processes. While the experiments, standing alone, cannot 
rule out other phenomena, they lead us to tentatively at- 
tribute a major role to the spin dynamics in the changing 
part of the attenuation. In order to get an insight on how 
this mechanism emerges, we resort to the further hints 
provided by numerical calculations on model systems. 

IV. NUMERICAL SOLUTIONS 

The ZME sequence can only change the sign of the sec- 
ular part of the homonuclear dipolar Hamiltonian. Then, 
all the non-secular interactions are uncontrolled processes 
that contribute to the attenuation of the PE. Here, we 
will consider their effects calculating the exact evolution 
of the LP during the defocusing and refocusing periods 
in a cyclopentadienil ring and in a complete molecule of 
ferrocene with a single 13 C. 

Let us consider that the initial state, with only the I\ 
spin polarized, evolves during a time tn with the Hamil- 
tonian 

H 1 =W + Hi+Hl. s .+His = H y + £, (5) 

where 

H J = -*YiB V u y E lV k ( 6 ) 

k 

and 

Wis=J>2J*S*. ( 7 ) 

k 

Here, W (described by Eq. (@)) and , are the secu- 
lar and non-secular parts of the /-spin dipolar interaction 
respectively. Tiis is the non-secular heteronuclear dipo- 
lar Hamiltonian. Taking into account ideal (7r/2) pulses 
the evolution continues for a time 4r/2 with — \2\H y . We 
explicitly neglect the heteronuclear interactions during 
this period, since it was verified that they do not pro- 
duce any visible attenuation in the observed time scale. 
In this way, the S'-spin appears only in (0) , and its effect 
can be taken into account by replacing His —> J2k ^fc-^fc- 
This substitution reduces to one half the computations 
involved. Figure 6 shows the PE amplitude as a function 
of tji. The Gaussian fitting of the ferrocene experimen- 
tal data (Fig. 4) is included for comparison. Even when 
these numerical results cannot reproduce the experimen- 
tal data, they show that the non inverted interactions 
produce a considerable destruction of the dipolar coher- 
ence which manifests in the decay of Mpe- The dipolar 
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coupling with the 13 C spin is the most relevant term to 
produce attenuation in the presence of spin dynamics, 
although in its absence, it does not produce a significant 
decay, as it is shown by the oscillating curve in Fig. 6 
where a single proton is present. The calculations for 
a ring of five dipolar coupled protons and for two rings 
show that the attenuation increases as more spins are in- 
cluded. This trend implies that, in the highly connected 
actual interaction network, a stronger attenuation, ap- 
proaching the experimental curve, should be expected. 
Conversely, the sequence of solutions going up in Fig. 6, 
can be associated to a progressive reduction of the dy- 
namics obtained through the application of the REPE 
sequence. In analogy with the inset in Fig. 4, it would 
reach the absence of decay when dynamics is totally hin- 
dered (n — > oo). The lesson to be drawn from this cal- 
culation is that a small non inverted term, although non 
relaxing by itself, if mounted on a complex spin dynamics, 
it is amplified and leads to strong dynamical irreversibil- 
ity. Thus any reduction of the complexity of the dipolar 
order (increase of n in our experiment) will decrease the 
amplification effect. While this slows down the polariza- 
tion echo decay, an asymptotic time scale is not reached. 
This is the case of the ferrocene system. 

Another kind of uncontrolled processes that could con- 
tribute to the PE attenuation are those producing an 
exponential irreversibility. As a model for these we can 
consider a fluctuating isotropic magnetic field. For sim- 
plicity, we assume it is only effective during the evolution 
in the rotating frame. Thus, the evolution of the density 
matrix during tn is given by the quantum master equa- 
tkmQ: 

v a— x.y,z j 

= -\[ h} ,p] -S(p-p(oo)). (8) 

Notice that in the absence of dynamics (TL V = 0) we 
would obtain (If) ~ exp(— t^/r^). As in the previous 
calculation, the period of refocusing is governed by the 
Hamiltonian — [2]H y . Fig. 7 shows the results for a single 
cyclopentadienil ring. Initially the decay follows an expo- 
nential law with a characteristic time r^. However, after 
a while, it switches to an exponential with a shorter char- 
acteristic time. The typical time for this change is con- 
trolled by the dipolar dynamics. This can be clearly seen 
by considering different time scales (h/d) for the dipolar- 
interaction (Fig. 7). However, in our small systems, the 
decay rate of the new exponential is not affected by this 
dynamics. Simulations on systems with different sizes 
show that this rate increases with the number of spins 
involved in the dynamics. Hence, we can expect that 
in unbounded systems, where the number of correlated 
spins increases progressively, the decay could not be rep- 
resented by any finite number of switches. In that sense, 



when the dipolar dynamics is strong as compared with 
the exponential process (r m fc < r^), the observed Gaus- 
sian tail can be thought of as a continuous switch between 
progressively stronger exponential decays. In the recip- 
rocal case (t0 < Tmb), the experimentally relevant decay 
occurs while the exponential process is dominant. 

The lesson learned from these numerical results, is 
that there is an emergent phenomenon in large many 
body systems through which spin dynamics favors the 
irreversible processes. This phenomenon is the quan- 
tum analogous to the one observed in a classical gas of 
rigid disks |l9| ]. The many body interaction creates an 
instability that amplifies the numerical rounding errors 
progressively limiting the efficiency of the time reversal. 
Similarly, for a local spin excitation in a lattice, evolv- 
ing according with a cellular automaton dynamics [pof , a 
single spin flip (not included in this dynamics) prevents 
the full reversal of the polarization to the original site 
producing its spreading. 

V. CONCLUSIONS 

We have studied the decay of the polarization echoes 
in different systems. The ZME sequence in its original 
version, or even better in its modified version providing a 
proper normalization, allows one to distinguish between 
a tentative classification which recognizes two kinds of 
systems. They are the dipolar dominated ones, which 
present a Gaussian component that can be reduced by 
the REPE sequence (like ferrocene and cobaltocene), and 
those with strong sources of relaxation where an exponen- 
tial decay is already manifested with the ZME sequence 
(cymantrene). 

The REPE sequence, designed to reduce the time scale 
of the spin dynamics by increasing the number n of 
pulses, is able to give further information on the first 
group by revealing the existence of underlying relaxation 
mechanisms for large enough n's as it is the case of cobal- 
tocene. Looking for an experimental test to rule out pos- 
sible experimental artifacts in these results, one could 
resort to alternative pulse sequences which also achieve 
the dynamics reversal with less demanding train pulses 
[fHHl. Further tests of the sensitivity of the REPE se- 
quence as a tool to characterize new materials are ex- 
periments in which a dipolar dominated system is pro- 
gressively doped with local centers of relaxation (e.g. a 
solid solution of cobaltocene in ferrocene). Besides, the 
REPE sequence can be readily adapted to be used under 
MAS by rotor synchronizing the pulses according to the 
prescriptions of ref . || . 

Because of the deep physical meaning of the Polariza- 
tion Echo, beyond its potential applications as a practi- 
cal tool, the ZME and REPE sequences open an excit- 
ing new field of study in fundamental physics. It is the 
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onset of irreversibility in interacting many body quan- 
tum systems. This is a particular area of the growing 
discipline of quantum chaos |^3[ whose theoretical devel- 
opment is still in its infancy. On one side the magni- 
tude of the polarization echo coincides with the most 
convincing signature of quantum chaos ^4|: the ex- 
ponentially fast decrease of the overlap between states 
evolved, from the same initial condition with slightly 
different Hamiltonians Ti and Ti + E. The connection 
is evident for a system of N interacting spins since we 
can write Mps = 2 N Tr (puPn+T,) — I7 which is a way 
to evaluate this overlap. This expression has close re- 
semblance to some definitions of entropy |^5| . Indepen- 
dently, -ln[MpE(tR)} is a measure of how the inefficiency 
in our dynamics reversal increases with time. When 
the polarization is a conserved magnitude within the ob- 
served time scale, MpE(tp) is a measure of the inverse 
of the volume ttt R occupied by the polarization at the 
PE time: -ln[MpE(tp)] oc ln[S! tj? ] which gives an "en- 
tropic" meaning for both the overlap and the polariza- 
tion echo. To further visualize this, one might assume 
that after a few pulses the shape of the distribution func- 
tion for the refocused polarization (e.g. P(x,t) = l/fl t 
for |x| < r2j/2 where t is an echo time) remains the same 
(i.e. the echo would represent a scale transformation). 
In this case — J P(x) ln[P(x)]dx grows proportionally to 
-]n[P{0,tR)} = -HMp E (tR)}. 

Thus, the Gaussian regime for the attenuation of the 
PE in our experiments seems to indicate an anomalous 
decrease (growth) of the overlap ("entropy") not pre- 
dicted within any of the present one body theories. Our 
experiments might be interpreted as an indication that 
infinite interacting many body quantum systems present 
an extreme intrinsic instability towards dynamical irre- 
versibility. The experiments, even with imperfections, 
are much ahead of the possibilities of the current theory 
and they have the leading word. 
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Figure 1: a) ZME sequence for refocusing the dipolar 
evolution of the X H spin polarization in the laboratory 
frame (see text). The solid and dashed lines show the 
main and secondary amplitude polarization pathways b) 
Alternative version of the ZME sequence where the re- 
focusing occurs in the rotating frame. The polarization 
maximum at £2 = t m provides an experimental point 
for data normalization. It allows a straightforward com- 
parison of the PE attenuation in different systems, c) 
Reduced Evolution Polarization Echo (REPE) sequence 
to control the dipolar spin dynamics. For a given time 
tft, setting t\ = tji/n and ti = [5] t\, the development 
of multiple spin order is reduced by increasing n. 

Figure 2: Normalized polarization echo amplitude 
as a function of t R for polycrystalline samples of cy- 
mantrene, ferrocene and cobaltocene, using the ZME se- 
quence. The lines are fittings to a Gaussian (ferrocene), 
exponential (cymantrene) and a product of both func- 
tions (cobaltocene). 

Figure 3: 13 C NMR spectrum of a ferrocene single 
crystal using the REPE sequence with: n = l, uiu/2tt = 
62 kHz, tc — 2 ms, t$ — 1 ms, td — 53 /is and t R — 8 /is. 
The two peaks indicate the presence of two magnetically 
non-equivalent sites. The on resonance peak corresponds 
to molecules with its fivefold molecular symmetry axis at 
approximately 20° with respect to the external magnetic 
field. 

Figure 4: Attenuation of the polarization echo in the 
ferrocene single crystal as a function of t R for n = 1 
(REPE sequence) . The line represents a Gaussian fitting 
with the characteristic time r^ b = (245±5) /is as the only 
free parameter. Inset: PE attenuation for progressively 
reduced dipolar dynamics, n = 1,2, 8, 16. No asymptotic 
regime is reached within the experimental time scale. 
The solid lines are Gaussian fittings yielding character- 
istic times r^ b = (335 ± 10) /is,r,^ b = (650 ± 20) /is and 
t^ = (1070±60) M s. 

Figure 5: Attenuation of the polarization echo in the 
cobaltocene sample as a function of Ir Tor n — 1,2,5,8,16 
with tt>i//27r = 56 kHz, tc — 85 /is, t$ — 150 /is, and 
td = 85 /ts. The experimental data show a clear crossover 
between a dominant Gaussian attenuation to an expo- 
nential one. The solid lines represent fittings of the whole 
set of data to equation (||) yielding a = (130 ± 10) /zs, 
b = (90 ± 10) f-is and r = (640 ± 20) /is. 

Figure 6: Numerical calculation of the attenuation 
of the polarization echoes considering the effects of the 
non-inverted interactions. The thick solid and dashed 
lines correspond to a cyclopentadicnil ring and to a fer- 
rocene molecule respectively. The Gaussian curve is the 
fitting of the ferrocene experimental data of Fig. 4. The 
attenuation is mainly caused by the destruction of the 
dipolar coherence induced by the heteronuclear dipolar 
coupling. Tiis does not produce relaxation in the absence 
of homonuclear dipolar interactions (oscillating curve). 



Figure 7: Attenuation of the polarization echoes for 
a cyclopentadienil ring considering an exponential irre- 
versible process (See Eq. (||)). The different curves rep- 
resent a progressive expansion of the dipolar time scale 
(H/d). Here d\ — 2610 Hz x 2-nh is the nearest neighbor 
dipolar coupling. The initial decay rate, 1/t^ — 500 Hz, 
switch to a stronger one after a while. The upper curve 
correspond to the exponential decay, exp(— tfl/r^), in the 
absence of dipolar dynamics. 
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